tj in the nuclear medium within a chiral unitary approach 



T. Inoue* and E. Oset 

Departamento de Fisica Teorica and IFIC, Centro Mixto Universidad de Valencia-CSIC 

CN ■ Institutos de Investigation de Paterna, Apdo. correos 22085, 46071, Valencia, Spain 

O 
O 
<N 

Abstract 



(N 



The self-energy of an rj meson in the nuclear medium is calculated in a chiral unitary 
approach. A coupled channel Bethe-Salpeter equation is solved to obtain the effective rj- 
N interaction in the medium. The base model reproduces well the free space tt-N elastic 
and inelastic scattering at the rj-N threshold or A* (1535) region. The Pauli blocking on the 
nucleons, binding potentials for the baryons and self-energies of the mesons are incorporated, 
including the rj self-energy in a self-consistent way. Our calculation predicts about — 54 — i29 
MeV for the optical potential at normal nuclear matter for an rj at threshold but also shows 
^ \ a strong energy dependence of the potential. 

Ph 1 Introduction 

i ' 

P+ 

The nuclear medium effects on meson properties are interesting and have been investigated ex- 
tensively. The data from relativistic heavy ion collision at CERN [|l], |2|, [| and at lower energies 
at Bevalac M, [5], may indicate either a lowering of the p meson mass or a large broadening of 
^ ■ its width. In the near future, experiments at GSI(HADES collaboration || [7]]) could clarify the 
situation by providing better statistics and mass resolution. The recent discovery of the deeply 
bound state of 7r~ in Pb || [| reveals an upward shift of the pion mass of about 20 MeV, and 
confirms the repulsive S-wave interaction of the 7r~ in nuclear matter already established from 



studies of the bulk pionic atom data [10], |TT, 12] . The low energy magnitudes of the r/- N interac- 



tion and the properties of the rj meson in the medium, are still open questions. The r\ mesic nuclei 
are expected to provide such informations and are searched at several facilities. For example, the 
7 Li(d, 3 He)6He, 12 C(d, 3 He^B and 27 Al(d, 3 He) 26 Mg reactions are investigated at GSI @. 

In this paper, we study the properties of the rj in the medium from a theoretical point of view. 
We evaluate the effective S-wave rj-N interaction in the medium in a chiral unitary approach. The 
self-energy of r] is obtained by summing the rjN interaction in the medium over the nucleons in 
the Fermi sea. The present work follows closely the method used in |14[] to evaluate the K nucleus 
optical potential in a self-consistent way. In the next section, we explain our approach. We show 
our input related with the nuclear medium in section 3. The results are given in section 4. Section 
5 is devoted to the conclusions. 
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Figure 1: Diagrammatic representation of the rj self-energy from S-wave interaction with nucleons. 



2 Chiral unitary approach to the 77 self-energy 

We want to obtain the self-energy of the 77 meson in nuclear matter at various densities p, as a 
function of the r\ energy k° and the momentum k in the nuclear matter frame. In this paper, we 
calculate it by means of 



U v (k°,k -p) = 2 



k p d 3 p n 
(2^)3 



T r/n (P°,P;p)x2 



where p n and kp are the momentum of the neutron and the Fermi momentum at density p 
respectively, and T vn (P°,P ; p) is the r/-neutron in-medium S-wave interaction, with the total 
4-momentum of the system (P°,P) in the nuclear matter frame, namely P° = k° + E n (p n ) and 
P = k + p n . Here, the isospin symmetry, T vp = T vn , is assumed and the amplitude is summed 
over nucleons in the Fermi sea as shown in Fig.[l]. We shall be concerned about the S-wave 77 
self-energy. At low energies of the 77 this part of the potential is largely dominant. 

We evaluate the in-medium amplitude T vn in a chiral unitary approach. For this purpose 
we follow the model for the free space n-N and coupled channels scattering of ref. |15], which 
reproduces the experimental data of ttN scattering up to energies above the iV*(1535) region. A 
similar chiral approach which covers a wider energy range in isospin 1/2, although with more free 
parameters, is also done in ref. [|16|1 . 

In ref. [O] the Bethe-Salpeter equation is considered with eight coupled channels including two 



7r7riV states, namely {n 



T]Tl, 



A°A, AT+IT, A°S C 



p, 7r + 7r n). The kernels for the 



meson-baryon two-body sector are taken from the lowest order chiral Lagrangians and improved 
by applying a form factor corresponding to a vector meson exchange in the t-channel. The kernels 
for ttN <-> ttttN transitions are determined so that they account for both the ttN elastic and 
ttN — > 7r7riV processes. That model reproduces well the the nN scattering amplitudes, especially 
in isospin 1/2, for the center of mass energy energies from threshold to 1600 MeV. It reproduce 
also the n~p — > rjn cross section at the region where the P-wave contribution is negligible. In 
this coupled channels approach, the model also provides the r]-n interaction in free space and 
generates dynamically the iV*(1535) resonance providing the width and branching ratios for its 
decay in good agreement with experiment, among them the r)N branching ratio which is quite large 
for that resonance. The agreement of the model with the different available data around A* (1535) 
resonance region and the adequate description of the properties of the resonance, in particular 
the strong coupling to the r/N state, give us confidence that the model is rather accurate to make 
predictions on the i]N — > 7]N interaction. Therefore we use this model as a base to take nuclear 
matter effects into account. 

Analogously to the free space case, T vn (P°, P ; p) is given by one matrix element of the matrix 
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Figure 2: Diagrammatic representation of the Bethe-Salpeter equation. 

equation in the space of coupled channels 

T(P°, P ■ p) = [l - V(^)G(P°, P ; p)} _1 VWa) (2) 

which is the solution of Bethe-Salpeter equation(Fig.^|), where the kernel matrix V have been 
factorized on shell as shown in ref.[|I7]] and ref.|18[, and G is a diagonal matrix of loop functions. 
The kernel has nothing to do with nuclear matter and does not depend on density and hence can 
be parameterized in terms of only the invariant energy yfs or s = (P ) 2 — P 2 . One can find the 
explicit form of it in ref . |H| . The effects of the nuclear medium are taken into account through the 
loop functions Gi(P°, P ; p), which describe the propagation of intermediate states in the medium. 
The meson baryon loop functions for the free space are given by the integral 

G (P° P) = I ^ Ml - 1 (3) 

lK ' ; J (2tt) 4 E t {P - q) P° - q° - E t {P - q) + it (q ) 2 - f - mf + le K> 

which is the zero density limit of the loop function in medium. The generalization to the case 
where the initial meson baryon system is not in the CM frame (P ^ 0), is a necessity in order 
to evaluate the meson baryon amplitude for arbitrary meson and baryon momenta in the frame 
of nuclear matter at rest, where the r] self-energy is evaluated (see eq.([TJ)). We regularize the 



integral by means of a cut off rather than by dimensional regularization as done in ref. ||15|| . The 
choice of a cut off is preferable when one performs the calculation in nuclear matter since Lorentz 
covariance is manifestly broken given the fact that one has a privileged frame of reference, the 
one where nuclear matter is at rest. The integration variable q is the momentum of the meson in 
the loop, but in order to obtain a Lorentz invariant quantity when p = and P ^ 0, the limits 
in eq.(||) must be implemented in such a way that the meson momentum in the rest frame of the 
original meson baryon state, q cm , should be smaller than the cut off taken. Another possibility is 
of course to make a boost and work in the meson baryon CM frame which is the option followed 
in ref. [|14| . Appropriate choices of the cut off and the subtraction constants(aj) give the equivalent 
loop functions as in dimensional regularization (see also ref.[pl|]). Evaluating the integral in CM 
frame of the meson baryon system, we find that the choice 

q™* = lGeV , a nN = 34MeV, a m = 14MeV, a KA = 39MeV, a xs = -22MeV (4) 

gives meson baryon loop functions equivalent to those in our previous paper. In order to find 
G(P°, P; p) which appears in eq.(0) we shall use the formulation of eq.@ and eq.(f|) in the free 
case but will take into account Pauli blocking in the intermediate nucleon states, plus the meson 
and baryon self-energy in the intermediate states. 

The Pauli blocking is one of the important nuclear matter effects. We incorporate this by 
replacing the free nucleon propagator as 

1 e(\P-q\-k F ) 



P° - q° - Ei(P - q) P° - q° - E[(P - q) 



(5) 
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with the unit step function 8(\P — q\ — hp), which prevents the scattering to intermediate nucleon 
states below the Fermi momentum. 

Another important nuclear matter effect is the dressing of hadrons. All the baryons and 
mesons in the intermediate loops interact with nucleons of the Fermi sea and their dispersion 
relations are changed. For the baryons, we incorporate this effect, within a mean-field approach, 
as a momentum-independent binding potential. Due to baryon number conservation, only the 
difference between the nucleon and hyperon potentials is relevant for our purpose and is introduced 
in the hyperon propagator in the kaon-hyperon loops. For the mesons, we incorporate the dressing 
effect by introducing the self-energy function Ui(q°,q ;p). The meson propagators are replaced, 
using the Lehmann representation, as 



1 



(q°) 2 ~ q- 



mt 



(q°) 



mt 



n/(g°,g;p) 



la, SlM '* p) 



(q 



0\2 



u 2 + ie 



(6) 
(7) 



where Si(q°, q ; p) is meson spectral function given by 



Si(q°,q ]p) = ~ 



Im[n,(g°,g»] 



n\(q°) 2 -q - mj - Ui(q° , q ; p)\ 2 



as a function of energy and momentum in the nuclear matter frame. 

Summarizing, the in-medium ttN loop function, for example, is calculated by 



G-kn(P° , P ', P) — CL-nN + 



(9) 



d A q 



with 



M 



N 



6{\P-q\-k F ) 



^2ivS n {uj,q-p) 



E N (P -q)P°-q°- E N (P -q)+ieJo (q ) 2 - u 2 + ie 
P-q 



1 



\P\ 



q 



P + q 



(10) 



from the Lorentz boost. 

The 2-loop function needed for the irirN channels is also modified in nuclear matter. We still 
take a zero constant for the real part according to our previous study in the free space [ 15| . On 
the other hand, the imaginary part is modified along the lines of G^n and we obtain 



lm[G n7TN (P°,P ;P)} 



(2^/(2^^-^ ^(\P- qi -q 2 \-k F 



x6 (P° - E N ) n JltvS^w, qi ; p)S w (P° -E N -u, q 2 ; p) 
with E N = E N (P — qi — q 2 ) when both the Pauli blocking and the dressed pion are considered. 



3 Input 

We use commonly accepted values of baryon binding potentials which are —70 p/p MeV for the 
nucleon and —30 p/p MeV for the hyperons. Hence, the difference +40 p/p MeV is added to 
the hyperon energy Ey in the hyperon propagator. 
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Figure 3: Diagrammatic representation of the P-wave(left) and S-wave (right) term of pion self- 
energy. The crossed diagrams are omitted in the figure. 



The effective K-N interaction in the medium is studied in ref.[20, |21], O and the resulting 
kaon self-energy is Uk(p) = 0.13 m 2 K p/p approximately. We use this real constant self-energy 
for kaons in the kaon-hyperon loops. 

The pion self-energy which we use in the 7rN and irnN loops, consists of a P-wave term and 
an S-wave term as shown in Fig|| The P-wave term includes N-h, A(1232)-h and iV*(1440)-h 
excitations, and is basically given by 

K(q°,q;p)= ( ^] f {U N {q\q ; p) + U A (q°,q ; p) + U N *(q°,q ; p)} (12) 



where f n7T jy = 1.02, and U^iq ,^ ;p), U A (q°,q ; p) and U]y*(q ,q ; p) are the Lindhard functions 
corresponding to the above excitations respectively. For example, [/jy. (g°, q ; p) is given by 

\ 2 



U N *(q°,q;p) = U 



fixNN* \ M N * 

Unn ) P \q\k F 

M N * 



1 / 1 2m Z + l 

Z+ 2 (1 - 2 )1 " — 



✓ + 5(1 



,/2n 



In 



z' + 1 
z' - 1 



\q\k F 
M N * 
\q\k F 



-Q° - ^ - ( M n* - M N ) + l -T N *(-q°,q[ 



(13) 



(14) 



(15) 



where we use f n NN* / fnNN = 0.477. This term is small compared to the A-h excitation term whose 
coupling is f^NA / fnNN = 2.13, but has a sizable contribution for energetic pions. The expression 
for U A is analogous to that in eq. ( |T3"1 - |15| ) and is given in ref. [^], the factor 3/2 in eq.(|l~3"D is replaced 
by 2/3 and the coupling, mass and width of the iV* are replaced by the corresponding magnitudes 



for the A. In the actual calculation we use the P-wave contribution eq.(12) improving with a form 
factor, a short range correlation and a recoil factor as 



ttNN 

777- 



NN 



g 2 F(g°,g) 2 5:^(g ,g;p) 



Ui(q°,q;p) 



EiUi(q°,q ;p) 



l ~ (^f) 9'T,Mq\q-p) 



(16) 
(17) 



1 



U N (q°,q ;p) 



2M N , 

l~2lL) 2 UN*( q ,<i;p) 



1 



i) U A (q°,q;p) 
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Figure 4: Self-energy (left), propagator (center) and spectral function(right) of a pion for 500 
MeV momentum, at normal nuclear matter density. The dotted and dashed dotted lines in the 
self-energy (left) correspond to the conventional ones including only the N-h and A-h excitations. 

respectively, where F(q°, q) is a monopole form factor with cutoff of 1.2 GeV and g' = 0.6. 
The S-wave term of the pion self-energy is taken in the so called "Tp" approximation, 

n^(g°, q;p)=T nN xp= |ir 1/2 + jjr 3/2 } x p (19) 

using the free space S-wave tt-N scattering amplitudes calculated in the present model. This term 
includes the contribution of iV*(1535)-h excitation through the isospin 1/2 amplitude Ty2, and is 
small compared to the P-wave one and becoming only sizable for very energetic pions with about 
700 MeV momentum. 

The total self-energy of pion is given by n^(g°, q ; p) +H^,(q°, q ; p) and in Fig.f| we show the self- 
energy (left), propagator (center) and spectral function(right) of the pion for 500 MeV momentum 
for normal nuclear matter density. In the self-energy(left) graph, the accompanying dotted and 
dashed dotted lines correspond to the case when we take only the N-h and A-h excitations. They 
show that the contribution of iV*(1440) is relatively small. 

One of the needed input in our evaluation is the self-energy of the rj which appears in the rjn 
loop. However, this is not known beforehand since this is what we want to calculate here. We 
perform for this purpose a self-consistent calculation. In a first step, the Tp approximation is used 
since we have previously calculated the free rjn T matrix. With this we obtain a new self-energy 
of the rj which is introduced in the second step, and so on until convergence is reached and one 
obtains an output r\ self-energy equal to the input one. 

4 Results 

Fig.|| shows the free space 7/-neutron scattering amplitude T m obtained in the present model. It 
is plotted as a function of the invariant energy y/s, which is 1487 MeV at the threshold. One sees 
that the amplitude changes drastically, even from attractive to repulsive, above the threshold. 
This is due to the coupling to the iV*(1535) resonance which is generated dynamically in the 
present model. Note that the threshold is about 50 MeV below the center of the resonance. We 
can obtain the rj self-energy as an explicit function of k° and k as independent variables. One 
of the interesting applications is the study of n bound states in the nucleus, where the rj will 
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Figure 5: The free space T^ n amplitude as a function of the invariant energy y/s. The arrow 
shows the rjn threshold. 



have small energy and momenta. For this purpose we can evaluate the rj self-energy for k° = 
and k = as a good approximation. Non-local effects from the consideration of the energy and 
momentum dependence of the self-energy give very small corrections for weakly bound states as 
demonstrated in ref.|23] for K~ atoms. 

In order to evaluate the rj self-energy in a first step, we substitute the free rjn amplitude in 
eq.(fj), setting k = 0, and obtain the rj self-energy for zero momentum shown in Fig|| top. This 
simplest calculation is nearly equal to the Tp approximation. The magnitude of the self-energy 
increases almost linearly in p and its shape reflects that of T m . 

Then we turn on the Pauli blocking in the way explained in section 2. The resulting r\ self- 
energy for zero momentum is shown in Fig|| middle. Comparing to the top one, we see the 
interesting effects of the Pauli blocking. First of all, the resonant shape is strongly enhanced. 
This can be understood as the reduction of the decay width of the resonance, since the Pauli 
blocking forbids the decay to nucleons with momentum smaller than kp. Second, the peak of the 
resonance, seen in the imaginary part of the amplitudes, or equivalently in the zero of the real 
part, is only shifted moderately to higher energies. This is in contrast with the K case, where the 
A(1405) resonance in the KN interaction, is about 80 MeV shifted upward [24, 14]. As pointed 



out in ref.[24], this is natural because a large part of the A^*(1535) resonance generated in the 



present model, comes from the K + Y1, K°TP and K°A components [^5|, |2T|, [Tj|, which have nothing 



to do with Pauli blocking. In ref. |24| it is argued that the Fermi motion compensates the decay 
width reduction due to the Pauli blocking. There also is some effect of this sort here, but we still 
find a net reduction of the A^*(1535) width due to the consideration of both effects, resulting in an 
enhanced strength of both the real and imaginary parts of the rj self-energy. A major difference in 



our model from ref. [^4j , at the present stage of the calculation, is the presence of the 3-body nnN 
channels. Also, our free space T vn (y/s) shown in Fig.[5], is already different from that in ref. ||24|| , 
where the peak of the resonance is seen at around 1590 MeV. 

Finally, we turn on the dressing of the hadrons in addition to Pauli blocking. This study is 
the main novelty of the present work. We use the input explained in the previous section. Besides 
this, we need the spectral function of the rj in order to calculate the rjn loop function. As stated 
in the previous section, we determine it in a self-consistent way by iterating the Bethe-Salpeter 
equation. This means that we need to study the finite momentum case even if we are interested 
only in the rj with zero momentum. We, hence, evaluate the rj self-energy also for finite momentum 
and, for example, Fig.|7| left shows a rough sketch of the rj spectral function obtained at normal 



7 



.15 



0.1 
0.05 

% o 

^ -0.05 

TO 

% -0.1 

t: 

-0.15 
-0 . 2 




-0.25 



Re 



0.3p 
0.5p 
0-VPo 
P±_ 



Im 



(a) Free 



-150 -100 -50 50 
*°-m, [MeV] 



100 -150 -100 -50 50 100 150 
k°-m n [MeV] 




Re 




Im 







(c)Full V/ 












7Po - 




Po ~ , 





. z j 

-150 -100 -50 50 100 -150 -100 -50 50 100 150 
fc°-m, [MeV] t°-m, [MeV] 



Figure 6: self-energy for zero momentum as a function of energy, for four different densities and 
in three different approximations. Top, (a) Free: the free space 77- N amplitude is used. Middle, 
(b) Pauli: the Pauli blocking is taken into account. Bottom, (c) Full: both the Pauli blocking and 
hadron dressing are taken into account. 



8 







Figure 7: Left: A rough sketch of the rj spectral function at normal nuclear matter density, where 
the peak value is not realistic because of the plotting resolution. Right: Self-energy of the 77 with 
400 MeV momentum at normal nuclear matter density as a function of energy. 

nuclear matter density as a function of k° and k. The results in what follows include the dressed 
77 accounting for the 77 self-energy as a function k° and k. 

The 77 self-energy for zero momentum obtained in this self-consistent approach is shown in 
Fig.|| bottom. One can see the effects of the hadron dressing in comparison to the middle figure 
of the panel. Apparently, the strength is spread wider. This is expected because the 77 spectrum 
is distributed as stated above. Even then, the resonant shape is still clearly visible. The center of 
the resonance is shifted upward, for example, about 25 MeV for normal nuclear matter density. 
This means that the sum of the repulsion on K and so on, surpasses the attraction on the pions. 
The difference of the hyperon binding energy and the nucleon one, works as a repulsion for this 
calculation. The dressed 77 also shifts it upward slightly because the rj with large momentum(larger 
than about 250 MeV) feels a weak repulsion, as shown in Fig.|7] right. Recall that the KY,, KA 
and rjn components are more important than the nN and tcttN components for the resonance in 
the present model. The 25 MeV shift for normal nuclear matter density is small(less than 2% of 
the iV*(1535) mass) and consistent with 77 photo-production data pq| , which together with the 



theoretical analysis of ref. J27j assuming no shift of the resonance, suggest that the in-medium mass 
of the resonance is almost the same as in free space. 

Fig.[8] left shows the spectral function of a zero momentum 77 corresponding to the final self- 
energy. We see still a narrow peak even for normal nuclear matter density and, hence, it makes 
sense to regard the 77 in the medium as a quasi-particle with a modified mass. It is interesting to 



note a second peak in the spectral function at higher energies, which, as noted in ref. j 24] 



comes 

from the coupling of the 77 to the iV*(1535)-h. We can observe in the figure that as the density 
increases the two peaks move away from each other. 

We define the effective mass m*(p) as the energy k° which satisfies 

(k o ) 2 -m 2 ri -Re[U v (k ,0;p)} = 0, (20) 

which is almost the position of the peak in the spectral function. The result is plotted in Fig|| 
right, which indicates that approximately one has 



n ~ 1 - 0.05— (21) 



m v p 



9 



-100 -50 50 



*°-m, [MeV] p/p„ 



Figure 8: Left: rj spectral density for the zero momentum. Right: Effective mass of the r\ as a 
function of density. 



in short. This moderate downward shift agree with the result of p4| |. While, a much stronger 



about 11%, downward shift for normal nuclear matter density, is reported in []2S| where the quark 
meson coupling model is used. It should be interesting to test these predictions from experimental 
data of 77 mesic nuclei. 

The optical potential U v (p) and the effective rj- neutron scattering length a* n (p) defined as 

U M = (22) 



O) = "TZ TFTZ ^"^"^ (23) 



M n U v (m v , 0,p) 
An M n + m v p 

are useful when we study the bound state of an rj in nuclei using the local density approximation. 
We plot them in Fig.|9|. One can see that the real and imaginary parts of the optical potential 
deviate as the density increases. The real part becomes larger than the linear Tp approximation 
while the imaginary part is nearly Tp. This behavior is better seen in the effective scattering 
length, where the real part increases and the imaginary part almost stays the same. The depth 
of the optical potential at normal nuclear matter density is —54 — i29 MeV. This value is not in 
contradiction with the 5% (28 MeV) mass shift of the rj, because the self-energy is evaluated at an 
rj energy equal to the rj mass, while the effective rj mass, k° in eq. (|20|) , is obtained using the rj self- 
energy calculated at same k° energy. It is interesting to compare these result with the ones obtained 
in the literature. In ref.|24]] the potential obtained was U v (p) ~ (—20 — i22) p/po MeV which 
provides an imaginary part similar to ours but the real part is about one half of the one obtained 



here. In ref.[]29| the potential obtained assuming that the mass of the iV*(1535) did not change in 
the medium, as we showed it was approximately the case here, was U v (p) = (—34 — z24) MeV at 
p = po, somewhere in between the two former results. 

In order to facilitate the task of making an accurate as possible prediction for the rj bound 
states in nuclei, we parameterize our results for the r] self-energy as a function of energy and 
density. Given the strong energy dependence of the U v (k°,0 ; p) seen in Fig|j] for k° — m v < 0, 
the consideration of this explicit energy dependence in the Klein-Gordon equation should be 
important. We can parameterize our results in the region —50 MeV < k° — m v < 0, as 

Re[U v (k°,0 ;p)} = a(p) + b(p)(k° - m v ) + c(p)(k° - m v ) 2 + d(p)(k° - m v f (24) 
Im[rF# ,0 ;p)} = e(p)+f(p)(k -m v )+g(p)(k - mri ) 2 + h(p)(k -m v ) 3 (25) 



10 



Figure 9: Left: rj optical potential as a function of density. The box and circle stand for the real 
and imaginary parts of the potential in the Tp approximation with the threshold T^ n corresponding 
to a scattering length of 0.26 + i0.24 fm [|TjJ. Right: rj-n effective scattering length as a function 
of density. 



with 



a(p) = 


-36200.3 p/po - 24166.6 p 2 /Po MeV 2 


(26) 


b(p) = 


-1060.05 p/po - 326.803 p 2 / p 2 MeV 


(27) 


c(p) = 


-13.2403 p/po - 0.154177 p 2 /Po 


(28) 


d(p) = 


-0.0701901 p/po + 0.0173533 p 2 / p 2 MeV -1 


(29) 


e( P ) = 


-43620.9 p/po + 11408.4 p 2 /p 2 MeV 2 


(30) 


f(p) = 


-1441.14 p/po + 511.247 p 2 /p 2 MeV 


(31) 


g(p) = 


-27.6865 p/po + 10.0433 p 2 /p 2 


(32) 


Hp) = 


-0.221282 p/po + 0.0840531 p 2 /p 2 MeV" 1 . 


(33) 



By using the Klein-Gordon equation, and substituting p — > p(r) in the spirit of the local density 
approximation, one can then obtain binding energies and widths of the rj state in different nuclei. 



5 Conclusion 

In the present paper we have used a chiral unitary approach successfully applied to the study of 
the irN interaction and its coupled channels, in particular the r)N channel, in order to evaluate the 
rj self-energy in a nuclear medium. We have taken into account the standard many body effects 
like Fermi motion and the Pauli blocking in the intermediate N states. In addition we have also 
included the self-energy of the baryons and mesons in the intermediate states, including the rj self- 
energy in a self-consistent way. The results obtained are interesting. While qualitatively similar 
to other ones obtained before, we obtain however a deeper potential with a real part about twice 
as big as in former, more simplified, studies. In addition we show that the energy dependence of 
the rj self-energy is very pronounced below the rj threshold and it would be interesting to consider 
it in studies of r\ bound states in nuclei. For this purpose we have parameterized our results in 
an easy analytical form which can be used to solve the Klein-Gordon equation for rj bound states 
or to interpret some physical processes where rj states close to threshold play some role [Q. The 
stronger real part of our potential compared to previous ones and the fact that the imaginary part 
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of the potential decreases rapidly as the r\ energy decreases, open new hopes that r\ mesic states, 
relatively wide, but narrow enough compared to the binding energy, could exist and be identified 
in actual experiments. 
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